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Abstract. The origin of relativistic electrons in the Crab 
nebula which are producing the broad-band flat radio 
spectrum of this prototype plerion has proved difficult to 
understand. Here I show that these electrons can be nat- 
urally explained as a relic population of the pulsar wind 
electrons that have lost most of their energy in the ex- 
panding nebula because of intensive radiative and adia- 
batic cooling in the past. The observed radio spectrum 
suggests that the initial slowing-down time of the pulsar 
was Tsd < 30 yr, which implies that it has been born with 
a spin period Pq ~ 3 — 5 ms, several times shorter than 
presently believed. Consistency of these results with the 
current data and the historical records is discussed. 

Key words: ISM: supernova remnants - ISM: Individual: 
Crab Nebula - radio continuum: ISM - pulsars: general 



1. Introduction 

Flat radio spectra, with power-law indices < < 0.3 
(for 5*1, oc ) represent one of the key features of the 
plerions (Weiler & Panagia 1978), that distinguish these 
fiUed-center supernova remnants (SNRs) from a typical 
shell-type SNR with the mean ar ^ 0.5 (e.g. Green, 1991). 
Meanwhile, the origin of the electrons producing these un- 
usual radio spectra is not yet understood (e.g. Green and 
Scheuer 1992, Woltjer et al. 1997). 

Radio electrons appear particularly enigmatic for the 
Crab Nebula (see Kennel & Coroniti 1984 - hereafter 
KC84) where the synchrotron spectrum extends with — 
0.3 (Baars et al. 1977) from 10^ Hz to 10^^ Hz, and then 
it steepens to ctopt ~ 0.8 in the IR/optical band (Mars- 
den et al. 1984) This requires an energy distribution of 
the electrons N{'~f) oc 7""'' with ac — 1 + 2ar — 1.6 ex- 
tending over 3 decades of energy (Lorentz factor) 7 from 
~300 to 3 X 10^, which then steepens to etc — 2.6. For the 
mean magnetic field i?* ~ 0.3 mC and the known age 
of the Crab Nebula, born in 1054, this steepening seems 
to agree with the radiative break Aa = 0.5 (Kardashev 
1962) both in the magnitude and position of the break. 
One might thus suppose that all relativistic electrons in 



the nebula, including the ones with 7 below the radiative 
break energy 7,br ~ 3x 10^ (at present), called "radio elec- 
trons" , represent a single population of particles currently 
accelerated at the pulsar wind termination shock (Rees & 
Gunn 1974 - RG74) that would produce an injection spec- 
trum (5(7) with ainj = 1.6 at 7 < 3 X 10^, and a-mj — 2.3 at 
higher energies (in order to explain the X-ray fluxes with 
ttx — 1.1 — 1.2). However, the mean electron energy in such 
a spectrum 7 < 10^, which is much less than the energy of 
the wind electrons upstream of the shock, 7*^ ^ 10^ (e.g. 
Kundt & Krotscheck 1980 - KK80, KC84, Arons 1996 - 
Ar96), whereas the energy and particle number flux con- 
servation laws across the shock require that 7 ~ 7*w An 
idea of current acceleration of radio electrons in the main 
nebula beyond the wind termination shock appears also 
problematic, because radio observations do not show any 
signiflcant variation of the radio spectral index, imply- 
ing an absence of effective electron acceleration sites there 
(Kronberg & Bietenholz 1992, Bietenholz et al. 1997). 

The remaining alternative is to assume that radio elec- 
trons in the Crab nebula are relics of the history of its 
evolution (Shklovskii 1977, KK80, KC84). The model pro- 
posed below reconsiders the evolution of the electrons in- 
jected into the expanding nebula at the pulsar wind shock. 

2. Results 

Previous studies of expanding plerions (Pacini & Salvati 
1973 - PS73; Reynolds & Chevafier 1984 - RC84) have 
assumed an injection spectrum Q(7,t) oc 7~"'"j with 
Q^inj = l + 2ar starting from 7 > 1. Thus, the origin of such 
a flat spectrum was not considered. In agreement with the- 
oretical predictions (KK80, KC84, Ar96, Atoyan & Aha- 
ronian 1996), here I assume an injection spectrum with 
a profound deficit of electrons below some 7c, for exam- 
ple in the form Q(7,t) cx x^(l -|-x)~^~"'°^ exp(-7/7,„ax) 
where x — 7/7c- For ainj ~ 2.3, the mean energy of elec- 
trons in this spectrum 7 ~ (8 — 10)7c. This defines 7w of 
the pulsar driven wind of electromagnetic fields and par- 
ticles which powers the plerion (RG74). The electrons are 
injected into the nebula with a power — rj^Lsd where 
Lsd = ^0/(1 + t/Tsd)^ is the spin-down luminosity, and 
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k = (n + l)/(n — 1) in terms of the pulsar braking index 
n. A significant, if not the dominant (r^e — > 1, KC84), frac- 
tion of Lsd is injected in the electrons, so 7w(0 oc Lgd/N 
where N(t) is the production rate of e"*" — e~ pairs by the 
pulsar. Very generally, N can depend mainly on the mag- 
netic field -Bn of the neutron star and its spin frequency 
fl(t) = r2o(l +t/Tsd)~^^^^~^^ ■ Because the basic time de- 
pendent parameter is f2, I approximate 7w oc QP{t) with 
a parameter p. This results in 7^ oc (1 + t/Tsd)~^ withQ a 
model parameter s = p/{n — 1). 

Because of intensive synchrotron losses in the past, 
electrons were brought to energies much smaller than 7*br 
at present. For the energy losses Pg = —{d^/dt)s = 67^ 
where b « 0.04(5/1 G)^ yr""'^, and magnetic field in the 
expanding nebula declining in time a,s B = B^{t^/t)"^, 
typically with m ~ 1.3-2 (RC84), the equation 7/Ps = t 
gives the radiative brake energy 



7br(t) = 3 X 10^(5,70.3 mG)-2 (t/t,)2 



(1) 



All electrons injected into the nebula at times t with en- 
ergies higher than 7br(i) are rapidly cooled down to this 
energy. Electrons of smaller energies are affected only by 
adiabatic energy losses which do not change a power-law 
shape of A^(7) once established. 

The origin of the flat distribution of radio electrons 
in plerions is readily understood if we neglect, for a mo- 
ment, the adiabatic energy losses. Then Eq.(l) predicts 
that all electrons that enter the nebula between times t 
and t + At with energies above 7bi (0 will be soon found 
in the energy region between 7 « 7br and 7 -I- A7, where 
A7 oc t^"^^^At. Since injection occurs mostly at high en- 
ergies, a (5-functional approximation Q{j, t) — ^¥^(7— 7w), 
where N = iVo(l-|-t/Tsd)'*~'^, can be used to find the num- 
ber of electrons AiV = NAt accumulated in the interval 
A7. Eq.(l) connects the injection time of electrons with 
their energy after radiative cooling as t oc 7i/(2m-i)_ Then 
the distribution AA^/A7 N{'-f) of electrons that have 
been injected at times t <^ Tsd acquires a power-law slope 
oc 7~"" with ao = 1 — l/(2m — 1). At time t r^d the 
injection rate N{t) starts to decline, resulting in an in- 
crease of the slope to ae ~ 1. This can explain the spectra 
of plerions like 3C 58 or G21.5-0.9 where a^- ^ 0, and 
suggests that such plerions are powered by pulsars with 
Tsd not much less than their age. However, radio spectra 
with tti- ~ (0.2 — 0.3) are produced by the electrons in- 
jected at times t ^ Tsd, when N{t) rapidly declines, and 
ae 1 + {k — s — l)/{2m — 1). This slope extends up 
to the current radiative break energy 7,bn if the energy 
of electrons in the wind presently 7,^ > 7,br, as in Crab 
nebula, and to a smaller energy otherwise. 

This qualitative analysis demonstrates the role of the 
radiative losses for the formation of flat spectra of ra- 

^ Note that this relation predicts 7„ ~ const for times 
t <SC Tsd. This seems to be correct also generally, allowing a 
reasonable variation (e.g. a gradual decline) of B-nit), although 
this would somewhat change the relation s = p/{n — 1). 



dio electrons in plerions. The results do not significantly 
change after we take into account also the adiabatic energy 
losses of electrons, Padb = avj/R. In the approximation 
of a homogeneous spherical source of a radius R which 
expands with a velocity v, the parameter a = 1 (Karda- 
shev 1962). For an arbitrary a, the energy distribution of 
electrons reads (Atoyan & Aharonian 1999): 



R 



(2) 



where R\ = R{ti), and Ff^ = r{j,t,ti) is the particle 
energy at an instant ti < t: 



r(7,i,ti)=7( ;|- 



1 



■(3) 



Expansion of the Crab nebula is approximated as v — 
f, (i/t,)* with acceleration parameter S w 0.1, as observed 
(Trimble 1968). In the region 7sd ^ 7 < 7*br, where 
7sd ~ 7*br (Tsd/i*)^™^^^""''"'' is the contemporary energy 
of the electrons injected at time t = Tsd, the 5-functional 
approximation for (5(7) results in N{'-/,t^,) oc 7"""= with 



Qfo = H- (fc - s - l)/(2m -1 + a + aS) 



(4) 



which confirms the result derived qualitatively above for 
a — 0. The spectral index a-^ = {ac — l)/2 then implies 



k -1- 2ar (2m -1 + a + aS) 



(5) 



The energy distribution of radio electrons in the 
Crab nebula should extend with the same index 
over 3 decades of energy below 7,br- This requires that 
10~'^7*br 3> 7sd, therefore the radio spectrum of the Crab 
nebula would suggest r^d <C l0-3/(2m-i+a+a5)^^ _ ^qq^^ 

(!)■ 

Numerical calculations confirm this very unexpected 
and far reaching conclusion. It may seem at the first glance 
to Fig. 1 that all curves explain equally well the measured 
synchrotron fluxes from radio to optical, and even fur- 
ther to X-ray and 7-ray regions (where the fluxes can be 
sensitive to spatially inhomogeneous distribution of the 
nebular magnetic field). However, the radio fluxes of the 
Crab nebula are known, and must be explained, with high 
accuracy < 10%. A closer look to the fluxes S'y shown 
on the inset in Fig.l reveals that an explanation of the 
data at low frequencies is acceptable only if Tsd ~ 30 yr or 
smaller. This implies that the spin frequency of the Crab 
pulsar initially was (1 -I- i*/Tsd)'^^*-"~^'' ~ 10 times higher 
than presently. For the spin period P = 27r/f2 this corre- 
sponds to the initial Pq « 3 ms. It could be a bit larger, 
Pq < 5 ms, for the canonical value of the pulsar braking 
index 71 = 3 (the radio data then require Tsd < 20 yr). In 
any case, the derived initial spin period of the Crab pulsar 
is much shorter than Pq ~ 19 ms predicted previously (see 
Manchester & Taylor 1977). Obviously, the reason for such 
puzzling discrepancy is to be understood prior to consider 
the model as realistic. 
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Fig. 1. The radiation fluxes (top) and the distribution of 
electrons (bottom) calculated for different Tsd- The pa- 
rameter s has been adjusted to give the best possible fit 
to radio data: s — —0.7 for Tgd — 20 yr (dashed); s = —1.1 
for Tsd — 40 yr (solid); s = — 2 for Tsd — 100 yr (dot- 
dashed) . The dotted lines are calculated in the formal case 
of Tsd = 1 yr assuming s = —0.3, in agreement with Eq.(5). 
Other model parameters are fixed at m = 1.4, a = 1, 
k = 2.3, S — 0.1. The fit at the optical and higher fre- 
quencies is reached assuming aiaj = 2.3, 7max = 2.2 x 10^, 
and 7c = 1.9 X 10^ (the principal parameter actually is the 
mean 7 « 1.6 x 10^). Compilation of data is from Atoyan 
& Aharonian (1996) and references therein. The size of the 
diamonds in the inset corresponds to a typical accuracy 
±6% of the fluxes reported by Baars et al (1977). Note 
that the negative values of s correspond to a gradual in- 
crease of 7w in time. Qualitatively, this can be explained 
as a result of very high efficiency of the cascade (Stur- 
rock 1971) of e+ — e~ pairs and 7-rays in the pulsar mag- 
netosphere. In the past, when the pulsar rotated faster, 
this cascade was developing in a more compact {oc c/fl ) 
and tense magnetosphere, where electrons were effectively 
produced in the thus much higher magnetic fields. Fast 
increase of N would then result in a smaller 7,^. 

3. Discussion 

This puzzle is solved if we remember that Pq has been 
previously derived from the relation 

t*=T*[l-(f^,/Oo)"-i], (6) 

with T* = f2*/(n — l)f2», which is found after integration 
of the torque equation O = —K f2" assuming that n and 
K do not depend on time. In that case the braking index 
n can be expressed through the spin frequency fi* and its 
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Fig. 2. The rotation period P{t) of the pulsar result- 
ing from numerical integration of the equation dX/dt = 
—KX^ for the normalized frequency X = assum- 
ing K(f) = Kq[\ + ^exp(— t/tK)] with different values of 
parameters A and Ik, and for different braking indices n: 
A = Q and n = 2.5 (dot-dashed); A = 4.1, Iyl = 170 yr 
and n = 2.25 (dashed); A = 6.2, = 155 yr and n = 2.25 
(solid). The 3-dot-dashed line is calculated for A — 4.1, 
tx = 170 yr, but assuming a time-dependent n{t) finearly 
decreasing from the initial Ug = 2.6 to = 2.25. In all 
cases the current value of the parameter j7*f2*/(f2*)^ = 
2.5 is provided. 

time derivatives at present as n = f2*f2*/(f2*)^. For the 
Crab pulsar the measurements result in n = 2.5 (Groth 
1975), so for ~ 900 yr and r* = 1640 yr Eq.(6) predicts 
Pq w 19 ms, which then implies also that Tsd — 700 yr. 

Reliability of these predictions, however, is question- 
able if we keep in mind that Eq.(6) is only an approximate 
formula for the pulsar age, and as such, it tells us that the 
model of a rotating magnetic dipole (Pacini 1967, Gunn 
& Ostriker 1969) does provide a good agreement with re- 
ality because > 0.6t* for any Qq. However, attempts 
to treat this formula as a strict equation appropriate for 
derivation of fig would overestimate its real accuracy, and 
can be largely misleading because it becomes practically 
independent of fio already at times when fi, < r2o/2. Ac- 
tually, the equation (6) has to push ^lo close to the current 
ri* in order to be able to formally compensate the differ- 
ence between r, and the known age i,. Instead of that, a 
mathematically correct approach is to attribute this dif- 
ference to the correction terms due in the approximate 
formula (6). For example, assuming that n = const hut 
K = K{t), we would derive = t*[1 — (rJ*/f2o)"~^] x C* 
where C* = t^K{t^)/ J^" K{t)dt. For K{t) gradually de- 
clining in time, e.g. because of decreasing magnetic field 
B^(t) of the neutron star [K oc S^), a correction term 
C* > 0.6 would readily allow any ^ Some formal 
examples are shown in Fig. 2. Thus, the knowledge of 17,, 
17, and 17, at present may not be sufficient for determina- 
tion of Oq significantly larger than 17*. Information about 
the initial parameters of pulsars can be found in the radio 
spectra of plerions which they create. 
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The rcsiilts found for the Crab nebula allow important 
conclusions about its history. The initial rotation energy 
of the neutron star, and hence also of the explosion that 
produced the star, reach the level > 10^"'^ erg. This con- 
firms the earlier suggestion of Chevalier (1977) that Crab 
Nebula was born in a normal Type II supernova event. 

Meanwhile, the total energy residing in the Crab neb- 
ula presently is only ~ 4 x 10^^ erg (e.g. Davidson & Fesen 
1985). Where is then the initial 10^^ ergs of the pulsar? 
The answer becomes clear if we take into account that 
this energy has been injected into the nebula mostly by 
electrons (r^e » 1, KC84) with a low-energy 'cutoff' (!) and 
the mean energy 7 ~ (t) that at times t -^t^, was much 
higher than 7br(i) down to which they have radiatively 
cooled on timescale At < t. Because of that, basically all 
of the injected energy has been radiated away, and only a 
tiny fraction £{t) ~ 7br/7w -C 1 of the pulsar spin-down 
power Lsd{t) was depositing in the nebula. The total en- 
ergy accumulated in the nebula by time t would be then 
W{t) oc ^^'"(l -I- t/TsdY~'^. Remarkably, because 7*w is 
still larger than 7*br, this energy might be increasing even 
currently if 2m — fc + s > 0. In particular, in the case 
of m = meq ~ 1.4 predicted by the model for the mag- 
netic field evolving close to equipartition with the elec- 
trons, 2m — fc + s>Oifs> —0.5. Numerical calculations 
confirm this unexpected result. The increase of W(t) at 
times t » Tsd (to be compared with RC84) could ex- 
plain the expansion of the nebula with a noticeable recent 
(Trimble 1968), or perhaps even current, acceleration. 

The principal suggestion of RG74 that the magnetic 
field energy in the Crab nebula is presently close to 
equipartition with We of the electrons not just acciden- 
tally, but because magnetic fields can effectively reener- 
gize the electrons when Wb exceeds Wc, becomes impor- 
tant for explanation of the observed small decline rate of 
radio fluxes, ~ (0.16 — 0.19) % per year (AUer & Reynolds 
1985, Vinyaikin & Razin 1979). For m = rrieq 1.4, the 
calculated decline rate 

- d In S^/dt « [m(l + a,) + 2aar{l + 5)]/U ■ (7) 

is in a reasonable agreement with observations only if 

a <C 1. Adiabatic reenergization of radio electrons in the 
regions of magnetic field compression could significantly 
compensate the expansion losses, effectively reducing the 
adiabatic parameter from a = 1 relevant for a simple ho- 
mogeneous model considered here, to much smaller val- 
ues in the case of a spatially inhomogeneous study. Note 
that the regions of magnetic field compression, like those 
around numerous dense optical filaments in the nebula, 
may be distinguished by an enhanced radio intensity, e.g. 
as diffuse radio filaments (Swinbank 1980), but impor- 
tantly, they will have practically the same ar — 0.27 as 
the rest of the nebula (Bietenholz et al. 1997), because the 
adiabatic energy gain (or loss) process does not change the 
power-law index of the electrons. 



For the Crab pulsar, a speed of rotation 10 times higher 
than presently implies a very high initial spin-down lumi- 
nosity Z/o ~ 10^^ erg/s. The pulsar thus was injecting into 
the young nebula an energy about 3 x 10^^ ergs per year, 
which was emitted mostly in the form of X-rays and 7-rays 
that could not be detected at that time. However, during 
the first years after the explosion, a significant part of this 
energy should be absorbed in the then opaque fast expand- 
ing shell of supernova ejecta between the nebula and an 
observer. This additional energy input into the thermally 
emitting shell can be the reason why the light curve of the 
supernova 1054 was declining much slower than that of a 
typical Type II event (Minkowski 1971). 

The high initial luminosity of the pulsar could have 
also a strong impact on the dynamics of the ejecta out- 
side the nebula. Interesting consequences for the possible 
progenitor of the Crab Nebula could be deduced from the 
high angular momentum of the pulsar, and hence of the 
progenitor star itself. Thus, the fast initial spin of the Crab 
pulsar as appears evident from the radio spectrum of the 
Crab nebula, suggests new approaches for understanding 
the origin and history of this prototype plerion. 
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